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Abstract 

We present near-infrared spectral measurements of Themis family asteroid (379) Huenna (D~98 km) and its 6 km satellite using 
00 3peX on the NASA IRTF The companion was farther than 1 .5" from the primary at the time of observations and was approximately 

5 magnitudes dimmer We describe a method for separating and extracting the signal of a companion asteroid when the signal is 
p I not entirely resolved from the primary. The spectrum of (379) Huenna has a broad, shallow feature near 1 //m and a low slope, 

characteristic of C-type asteroids. The secondary's spectrum is consistent with the taxonomic classification of C-complex or X- 
^ pomplex. The quality of the data was not sufficient to identify any subtle feature in the secondary's spectrum. 

^'■Keywords: ASTEROIDS, SPECTROSCOPY 
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1. Introduction 

Since the discovery of the small moon Dactyl around (243) Ida 
by the Galileo spacecraft more than a decade ago, many more sys- 
tems have been d isco vered by either indirect methods (lightcuves, 
k g. [Pravec et al.l [20o3) o r direct detection (direct imaging or radar 
echoes, e.g. 



Merlineet all 19991 : iMargot et al.]|2002l : lOstro et al.ll2006l: 



iMarchis et alJl2008h . The study of the physical and orbital character- 
istics of these systems is of high importance, as it provides density 
measurements and clues on fomatio n processes through models (e.g. 
Iburda et al.ll2004 IWalsh et al.ll2008h . Primary formation scenarios in- 
clude capture (from collision or close approach), and fission and mass- 
loss from either YORP s pin up for sm all binary asteroids (D<20 km, 
IPravec and H anis 2007; Walsh et al. 2008 ) or from oblique impacts 
■for large asteroids (D>100 km, Descamps and Marchisll2008l) . 

Spectroscopic measurements of each component of multiple sys- 
tems are fundamental to understand the relationship between small 
bodies and their companions. From high quality measurements, dif- 
ferences in slope and band depths provide information on composi- 
tion, grain sizes, or effects of space weathering on each surface and 
help constrain formation scenarios. Obtaining photometric colors or 
spectra of the secondaries, however, is rarely possible due to the small 
(sub arcsecond, often less than a few tenths of an arcsecond) separa- 
tion and relative faintness (often several magnitudes dimmer) of the 
satellites. Photometric colors of TNO binaries, which are typically of 
comparable size, have revealed their colors to b e the same and are in - 
distinguishable from the non-binary population jBenecchi et alj200^ . 



'Visiting Astronomer at the Infrared Telescope Facility, which is operated 
by the University of Hawaii under Cooperative Agreement no. NNX-08AE38A 
with the National Aeronautics and Space Administration, Science Mission Di- 
rectorate, Planetary Astronomy Program. 



The first spectroscopic measur ements of a resolv ed asteroid system, 
(22) Kalliope, are presented bv^ Laver et al.l ( l2009 t) using the Keck AO 
system and its integral field spectrograph OSIRIS. They find the spec- 
tra of the primary and secondary to be simi lar, suggesting a common 
origin. More recently iMarchis et al.l (l2009h presented a comparative 
spectroscopic analysis of (90) Antiope based on data collected using 
the 8m-VLT and its integral field spectrograph SINFONI. This binary 
asteroid is composed of two large (D~85 km) components. A compar- 
ison of Antiope's component spectra suggest that the components are 
identical in surface composition and most likely formed from the same 
material. 



In November 2009, main-belt asteroid (379) Huenna reached per- 
ihelion, and the system h ad a separation of over 1.5". The pri- 
mary is a large (D~98 km. Ma rchis et al.ll2008l : iTedesco et al.ll2002h 
C-type (Bus a nd Bip zel 2002) and a member of the Themis family 
(Zappala e t^i l99j) . Its secondary (S/2003 (379) 1) was discovered 
by Margoli tOO^ With an estimated diameter of 6 km, the sec- 
ondary is ~5 magnitudes fainter at optical wavelengths than the pri- 
mary, and it has a highly eccentric orb it (e ~ 0.22) with a period of 
87.60 ±0.026 days (Marchis et al.ll2008t) . This irregular orbit suggests 
mutual capture from disruption of a parent body or even an interloper 
dMarchis et alj|200^ . It is not likely that the secondary formed by 
fragmentation of Huenna due to the YORP spinup effect, b ecause the 
primary is too la rge for this process to be efficient (Rubi ncamll200(]| : 
IWalsh et"aDl2008ilV Here we present near-infrared spectral measure- 
ments of the primary and secondary and as well as a new spectral ex- 
traction method to separate them. 
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2. Observation and Data Reduction 



2.1. Spectroscopy 

We observed the (379) Huenna system on 2009 November 28 UT, 
and obtained near-infrared spectral measurements from 0.8 to 2.5 //m 
using S peX, the low- to medium-resolution near-IR spectrograph and 
imager jRavner et al.li200 3). on the 3-meter NASA IRTF located on 
Mauna Kea, Hawaii. Observations were performed remotely from 
the Ce ntre d' Observation a Distance en Astronomie a Meudon (CO- 
DAM, iBirlan et al.ll2004h i n Meudon, Fran ce. Using the Binary Orbit 
Fit (BOF) algor ithm from toescamps (2005) and the orbital parame- 
ters provided in lMarchis et al. (2008.) , the secondary was expected to 
be located about 1.88" W and 0.17" N with respect to the primary at 
the time of our observations. The configuration of the slit (shown in 
Fig. [T} was thus chosen as following. The slit was first centered on 
the primary (position "P" in Fig.[T) and 8 spectra of 120 seconds were 
taken, alternated between two diflrerent positions (usually noted as 'A' 
and 'B') on a 0.8"xl5" slit aligned in the north-south direction. We 
chose not to align the slit along the parallactic angle to facilitate the 
positioning of the slit on the secondary for these challenging observa- 
tions. Next, the slit was shifted west ~1.8" (15 pixels; position "S" in 
Fig. [TJ to maximize the flux from the bulge on the southwest portion 
of the guider image, the secondary, and minimize contamination from 
the light of the primary. We continued to guide on the primary while 
taking spectra of the secondary. Because of the long period (~87 day) 
of the secondary, its motion relative to the primary during our observa- 
tions was negligible. We spent two hours taking spectra of 120 seconds 
of the secondary, although we combine and present only the 1 1 spectra 
toward the end of the night, when the seeing was lowest and the signal 
of the sec ondary was det ectable. The solar-type standard star Landolt 
102-1081 ( lLandoltl[l983l) was observed to remove the solar contribu- 
tion from the reflected light, and an argon arc lamp spectrum was taken 
for wavelength calibration. 

The dome was closed for the first two hours of the observations due 
to high humidity. During the beginning of the observations the seeing 
was around 1" and humidity was high. The data we present here, how- 
ever, were taken toward the end of the night when the humidity had 
dropped to 10% and the seeing was down to about 0.5". 

Reduction was performed using a combination of routines within 
the Image Reduction and Analysis Facility (IRAF), provided b y the 
National Optical Astronomy Observatories (NOAO) (lTodvlll993l) . and 
Interactive Data Language (IDL). A software tool called Autospex 
was used to streamline reduction procedures to create a bad pixel map 
from flat field images, apply a flat field coiTection to all images and 
perform the sky subtraction between AB image pairs. The primary 
spectra were combined and the final spectrum was extracted with a 
large aperture centered on the peak of the spectrum using the apall 
procedure in IRAF. The spectral extraction procedure of the satellite is 
described in Section 3. A wavelength calibration was performed using 
an argon arc lamp spectrum. The spectra were then divided by that of 
the standard star to determine the relative reflectance and normalized 
to unity at 1.215 yum. 

2.2. Imaging 

On the same night (see Table [T]l, we also observed Huenna using 
the Adaptive Optics system of the Keck-II telescope and its NIRC2 
infrared camera (1024x1024 pixels with a pixel scale of 9.94 x 10"^ 
arcsec/pixel). Two Kp-band (bandwidth of 0.35 1 yum centered on 2. 124 
^m) images made of 4 coadded frames each with an integration time 
of 40 seconds were basic-processed using classical methods (sky sub- 
traction, flat-field correction and bad pixel removal). In the resulting 




Figure 1 : Observing configuration: P marks the position of the slit for observa- 
tions of the primary using the IRTF observations. S marks the position for the 
secondaiy which was chosen based on the expected position of the secondary. 
Top: Adaptive-optics image obtained with NIRC2 at the Keck telescope. Bot- 
tom : Image obtained at the IRTF with SpeX slit- viewer, displayed at the same 
scale. Both images are oriented with north up and east left. The plate scale of 
the IRTF spectrograph and IRTF guider (slit- viewer) are 0.15 and 0. 12"/pixel, 
respectively. 
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Table 1 


: Observations on 


28 November 2009 


Target 


Time (UT) 


Airmass 


Instrument 


(379) Huenna 
S/2003 (379) 1 

L102-1081 
(379) Huenna 


12:56-13:08 
13:13-13:56 
14:19-14:26 
13:46-13:52 


1.31-1.37 
1.40-1.75 
1.34-1.31 
1.6-1.75 


Spex, IRTF 
Spex, IRTF 
Spex, IRTF 
NIRC2, Keck-II 



frame shown in Fig. [T] (averaged observing time 13:49:05), obtained 
after shift-adding these two frames, the primary is not resolved and 
slightly saturated. The satellite is detected with a peak SNR~23 and 
no motion with respect to the primary was detected over the 6 minutes 
of observation. 

The position of the satellite as measured on the NIRC2 image 
(Fig.IB is 1.67" ±0.05" W and 0.62 " +0.05" S, sig nificantly ofi'set from 
the predicted position based on the lMarchis et al.i(2008) orbital model 
(1.88" W and 0.17" N). This highlights the need for a refined period 
and adjustment for precession effects due to the irregular shape of the 
primary. Additional astrometric measurements are necessary to bet- 
ter constrain the orbital period and precession. For the spectroscopic 
observations on the IRTF, because the slit was oriented N-S and the 
displacement of the satellite in the E-W direction with respect to its 
expected position was only 0.2" (compared to a 0.8" slit), the slit was 
well-positioned on the satellite (see Fig.[TJ. 
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Figure 2: A cut of the (379) Huenna spectrum along a wavelength-constant 
line (/t = 1.23 jum). Acquired data (observing time 13:58 UT) is plotted as 
a black solid line, and the corresponding best-fit curve is drawn as a dotted 
red line. Individual Gaussian functions for Huenna (primary) and its satellite 
(secondary) are drawn in bulleted and solid blue curves, respectively. Central 
positions (Cp, c,), peak amplitudes (fp, /,), and Gaussian standard deviation 
(o"), accordingly to notation in Eq[T] are labeled for clarity. Residuals of the fit 
are plotted on the bottom panel. 



3. Spectral Extraction of the Satellite 

Given the spatial resolution of SpeX, which is limited by the at- 
mospheric seeing, and the angular separation of (379) Huenna and its 
satellite at the time of our observations (~1.5"), the spectrum of the 
secondary was contaminated by diffuse light from the primary. We 
thus extracted both spectra (the companion and the diffuse light of the 
primary) at once, by adjusting a fit function at each wavelength com- 
posed of the sum of two Gaussian functions (one for each component) 
on top of a linear background: 

fit = , . exp — — — + . exp — 

\ 2tr2 ; V2^ \ 20-2 ) 

where x describes the pixels perpendicular to the wavelength disper- 
sion direction. The 7 free parameters are: a and b, the slope and ab- 
solute level of the background, fp and /j, the maximal flux of each 
component centered on c,, and c,, respectively, and a the Gaussian 
standard deviation. In Fig.|2]we present an example of the adjustment 
showing both components and their sum plotted with the data. 

We verified the accuracy and robustness of the extraction by exam- 
ining the residuals of the fit and the evolution of the free parameters 
according to wavelength shown in Fig. [3] First, the angular separation 
between the two components, the peak of the diffuse primary light and 
the secondary, (f,,-Cj) is nearly constant, with a value of 1 .04 ± 0.25" , 
the uncertainty being mostly due to the uncertainty of the position of 
the secondary (c,, see Fig.[3j- This separation is the order of magni- 
tude expected (0.62" along the North-South direction, see Fig.[T). 

Second, the Gaussian standard deviation (cr) follows the expected 
trend due to the seeing power-law dependence with wavelength. In- 
deed, the seeing is proportional to /i"''^, which means that the amount 
of diffused light from Huenna entering the slit is not the same at short 
{e.g. J) and long {e.g. K) wavelengths. The same effect occurs with the 



slit-dispersed light: in the direction perpendicular to the wavelength 
dispersion, the spectrum is broader at shorter wavelengths and nar- 
rower at longer wavelengths. The combined effect is a /i^^/s depen- 
dence of the standard deviation (cr). We plot in Fig. [3] the theoretical 
wavelength dependence function of the Gaussian standard deviation, 
labeled "seeing", for a seeing value of 0.5" at 0.5 yum. 

We set the standard deviation of the satellite Gaussian function 
equal to Huenna's to simplify the fitting algorithm. However, from 
the discussion above, we can expect a wavelength dependence in /t"''' 
for the satellite (which is entirely covered by SpeX slit). This implies 
•we slightly underestimate the size of the satellite Gaussian function 
toward longer wavelength, and thus the satellite flux, and therefore its 
spectral slope. 

Third, the peak-to-peak amplitude measured for H uenna and its 
satell ite corresponds to the value of Am a; 6 reported bv lMarchis et al.l 
(2008): the ratio fjfp measured for the spectra is 0.4 ±0.1, and the 
amount of diffused light from Huenna entering the slit is about 1 % of 
the overall peak amplitude of Huenna (/«), calculated using a Gaus- 
sian profile for Huenna due to ~0.5-0.6" seeing, conditions similar to 
what we had at the end of the night. Note that is Huenna's overall 
maximum, not the local peak we measure for dispersed light that we 
call fp. 

Finally, the background level {b), mimics the spectral behavior of 
Huenna, whose diffused light dominated the light entering the slit, as 
expected. The amount of diffused light entering the slit 1" from Cp (a 
separation corresponding to 2cr) is about 10"' xfn, corresponding to 
K fp/250 (under the same seeing assumptions). 

We also test our extraction method by comparing the high quality 
spectrum of the primary while centered on the primary extracted by 
both methods. In the IRAF reduction, the center peak is identified, an 
aperture width is chosen, and all the flux within that region is summed. 
In our method we fit Gaussian functions to the flux. An additional error 
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is introduced in our extraction method because the fit functions can 
sUghtly over or underestimate the flux. The largest error is introduced 
before 1 /im and past 2.3 yum as is shown in the top panel of Fig. (3] The 
spectra are generally consistent, although the slope of the spectrum 
from the new extraction method is lower (0.70±0.09 vs. 1.04±0.05) 
and we do not see the 1 micron band as clearly because of the poorer 
fit of the function at shorter wavelengths. This suggests that the errors 
in the slope from the Gaussian extraction are greater than the formal 
eiTors listed. The spectra of the diffuse primary and the secondary have 
a low signal-to-noise ratio and thus high errors due to the low flux and 
the limited time when the seeing was sufficiently low. 



4. Results 

The spectra of Huenna and its companion are presented in Figure|4] 
Slopes are calculated between 0.9 and 2.4 fim using a chi-square min- 
imizing linear model (using the IDL routine linf it) and are listed in 
Table [2] with their formal errors for various wavelength ranges. The 
telluric regions from 1.3 to 1.4yum and from 1.8 to 2.0 pm were ex- 
cluded from the calculation. A spectrum of the primary published by 
IClarketal.l(l20rO ) is taxonomically in agreement with our spectrum, 
also a C-type, however the band minima positions differ (1.020/jm 
versus 1.045 fim for our data). The slope of our spectrum is also lower 
than the previously published data, although they agree within the er- 
rors, 1.04±0.55 and 1.94±0.53 %/10^ A. An additional eiTor of 0.5 
%/10' A is added to these slopes because in this case we are compar- 
ing data from different nights using different standard stars. 

The spectrum of the secondary is consistent with a C-complex or X- 
complex spectrum. The C- and X-complexes are broad umbrella terms 
for B-, C-, Cb-, Cg-, Cgh-, Ch-types and X-, Xc-, Xe-, and Xk- types, 
respectively. Without accompanying visible data, C-complex and X- 
complex spectra are indistinguishable except in the case of the pres - 
ence of a broad 1 fim feature indicating C-type jPeMeo et al]|2009l) . 
As can be seen in Table|2]the slopes of the spectra of the secondary and 
the diffuse light of the primary are greater than the slope of the primary 
spectrum even considering the large error. Particularly, the slopes are 
much greater in the shorter wavelength region, while the slopes are 
all relatively consistent in the higher wavelength region. Because this 
slope effect is seen for both spectra, we suspect it is not due to any com- 
positional differences but rather to the observational circumstances. 

Differences in slope could be caused by atmospheric difl^erential re- 
fraction which varies as a function of wavelength and is stronger at 
shorter wavelengths (Filiroenko 1982). For our observations, the slit 
was oriented in the north-south direction rather than along the parallac- 
tic angle. We chose this orientation to facilitate the positioning of the 
slit on the secondary while reducing the contamination of diffuse light 
from the primary. However, this slit orientation could introduce slope 
effects. Airmass differences at the times of the observations are likely 
the biggest contributor to slope effects. The airmass was 1.3 when 
centered on the primary and 1.4-1.7 when centered on the secondary. 
Because the airmass of the standard star was 1.34 at the time of obser- 
vation, effects of atmospheric refraction were not entirely removed for 
the observations when centered on the secondary. 

Considering the effects of airmass and parallactic angle that intro- 
duce an additional slope affect to the data, we are unable to be certain 
whether the slope of the secondary is consistent with the primary. This 
allows several possible formation scenarios. The satellite could be a 
fragment from a former parent body that was disrupted to form Huenna 
and the secondary. The satellite could be a member of the Themis fam- 
ily that was captured into an eccentric orbit around Huenna. Although 




0.5 1.0 1.5 2.0 2.5 
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Figure 3: Parameters of the spectral extraction plotted as function of the wave- 
length. The gray areas represent the wavelength ranges affected by telluric ab- 
sorptions, and we labeled the atmospheric windows with their usual designa- 
tion (visible, J, H, and K). From bottom to top. Sub-panel 1: position of the 
two Gaussian centers, Huenna (f and its satellite (c,). Sub-panel 2: the Gaus- 
sian standard deviation (cr) compared to the atmospheric theoretical seeing (see 
text). Sub-panel 3: the background level (b) compared to a fraction of Huenna 
flux (fp). Sub-panel 4: the residual of the fit, displayed as percent of the peak 
value of the spectra (/,, and /, for Huenna and its satellite respectively). 
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Table 2: Slope Measurements 



Object' 


Slope 0.9-2.4 fim 

[%/ilO' A)] 


Slope 0.9-1.3 yum 
[%/(103 A)] 


Slope 1. 4-2.4 yum 
[%/(103 A)] 


(379) IRAF 


1.04±0.05 


0.70+0.02 


1.29+0.02 


(379) Gauss 


0.70+0.09 


-1.20±0.05 


0.86+0.05 


(379) Diffuse 


2.34+0.59 


2.53+0.38 


1.45+0.38 


S/2003 (379) 1 


2.51 + 1.17 


3.41+0.77 


1.20+0.77 




*5 0.0 E i 

0.5 1.0 1.5 2.0 2.5 

Wavelength (fjm) 



Figure 4: Top section: From top to bottom: The spectrum of the primary with 
the slit centered on the primary reduced using IRAF; the same data reduced 
using the new Gaussian method presented here; the spectrum of the primary 
with the slit centered on the secondary measuring the diffused light; and the 
spectrum of the secondary. The spectra are plotted with one sigma errors and 
are normalized to unity at 1.215 /jm and are shifted by -0.5, +1.5, +2. 5, and +2.8, 
respectively. Bottom section: The ratio of the secondary and diffuse primary 
spectra both with the slit centered on the secondary. 



more unlikely, it could also be a captured interloper of similar taxo- 
nomic type. We are confident, however, that the satellite is not of a 
significantly different taxonomic type. 

We have presented the third resolved spectrum of a moonlet satel- 
lite of a m ultiple asteroid system, the first being (22) Kalliope from 
iLaver et al . (2009) and the second, (90) Antiope (Marchis et al. 2009). 
Our innovative extraction method allows us to derive the spectra of the 
satellite (with a magnitude difference of 5 and distance of 1.5") and 
the primary, and this method could be used for similar future projects. 
These results show the capability of a medium-size telescope to aid 
in the measurement of multiple system asteroids with large separa- 
tions. Further observations with larger telescopes, particularly those 
equipped with adaptive optics, can be used to perform the similar mea- 
surements for multiple asteroid systems which are not easily separable 
(angular distance <1"), and will be useful to constrain formation sce- 
narios. 
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